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ABSTRACT 


The real problem in acceptance sampling is to design or choose 
sampling procedures that accurately represent the consequences of 
accepting or rejecting the lot- Host of the acceptance sampling 
procedures in use today are based on slctltsitccil (probability) 
measures, with the decision maker assigning desirable 
probabilities to the events of accepting the lot at LTPD and of 
rejecting the lot at AQL- Sequential sampling, a valued approach 
when the cost of inspection is high, additionally aims at 
minimizing the expected number of items sampled. 
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loss of utility <of which monetary loss is a special case) to the 
decision maker in the sense, so that the plans 
obtained minimize the total expected risk of accepting/rejecting 
the lot, making full use of the existing (prior) knowledge of the 
unknown lot quality. The decision is guided by on one hand the 
loss of accepting or rejecting the lot, which directionally 
decreases with the number of items sampled (as more information 
of the lot quality is available) and on the other hand the cost of 
sampling, which increases with this number. 



As r isk~optimized sequential sampling plans cannot easily be 
obtained analytically, iterative techniques and algorithms are 
presented to obtain these plans- Numerical examples of such plans 
are given, and the results of sensitivity analysis are presented- 
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NOTATIONS 


a . 

1 

a . 

1 

ac 

ASN 

c 

C€«L I ( X ) 

cn 

6 

/ I aar ( x ) 
k 


LTPD 

n 


( n , r ) 


(n jr r ) 

iTi m 

fM 


N 




\P 


Prior probability of occurrence of - 
Posterior probability of occurrence of p^. 

Subscript designating the decisions "accept lot". 

Average Sampling Number of the plan. 

Cost of sampling one item. 

The smallest integer greater than or equal to x. 

Subscript designating the decisions "continue sampling". 
Loss due to imperfect information. 

The largest integer lesser than or equal to x. 
the number of components of the discrete k—point 
prior distribution « 

Lot Tolerance Percent Defective. 

Number of items sampled. 

Cumulative number of items sampled sequentially to reach the 
meeting point. 

Coordinates of a point in the sample space. 

Coordinates of the meeting point in the sample space. 

I (n ) . 
m 

The largest possible sample size (the effective meeting 
point) for a plan. 

The probability of an inspection station making a Type I 
error (rejection of a non— def e c t i ve item). 

The sampling path, a locus of (n,r) points as sampling 


proceeds . 



p 


The true unknown lot quality, the measure of the fraction of 


P 



r 


r 

fii 


defective items in a lot» 

The apparent lot quality as observed (will in general be 
different from p due to imperfect inspection). 

The i— th probable value of the distribution of the lot 
quality p. 

Probability of finding a defective item at the current 
stage of sampling. 

The "rejection count" — number of items found defective in 
sequentially sampling n items. 

Cumulative number of defective items found when the 
meeting point is reached. 


r (n ) 
ac 

r . (n ) 
rj 


rj 


R (n , r ) 
xy 


The rejection count at the acceptance boundary at n- 
The rejection count at the rejection boundary at n. 


Subscript designating the decisions "reject lot". 

The risk of making decision xy r j , cnj- j at a point 


R^(n , r ) 


(n,r) on the sample space. 

The optimum Bayes* risk at a point (n^r) on the sample 


e 

w ( p> 
ac 

y^(p) 

W^j(p) 

w^(p) 


space . 

The probability of an inspection station making a Type II 

error (acception of a defective item). 

Loss caused by accepting the lot at lot quality p. 

(See W (p) . > 
ac ^ 

Loss caused by rejecting the lot at lot quality p. 

(See W . ( p) . ) 
rj 



GLOSSARY^ 


AQL (Acceptable Quality Level) 

(ASQC Standard) 

The AQL is the maximum percent defective that for the purposes of 
sampling inspection can be considered satisfactory as a process 
average - 


ASN 

Average Sampling Number- The average number of items sampled in 
using a sampling plan- 

BAYES^ RISK (synonym for EXPECTED RISK) 

DECISION BOUNDARY 

The boundary on the ^airtpl^ spaces which separates the "continue 
sampling" region from one of the d,&ctsiori regions - 

EXPECTED RISK (of the sampling plan) 

The expected consequence of operating the sampling plan in terms of 
lass of the (expected) taT-nitriCLl decision- Called Bayes® Risk in 
literature- (See LOSS,) 

LOSS (caused by a decision) 

A quantified measure of all unfavorable consequences of a decision- 


] 



LOT QUALITY 


The fraction of defective items in a lot. 

LTPD (Lot Tolerance Percent Defective) 

The percent defective at which the consumer's risk is specified- 
The LTPD is considered to be a non-acceptable process average- 

MEETING poihnr 

The point on the spcLC:& where two d^cisic^n baurtd€Lr-t&i£ meet. 

(Wetherill £192). 

OC CURVE (Operating Characteristic Curve) 

The relationship between the probability of acceptance of the lot 
(based on a sampling plan) and the lat a^ucLlity, 

OSP (Optimal Sampling Plan) 

A sampling plan obtained by minimizing the &>cp&ct&d (or Bciy&s* ) 
r* t S'h « 

PRIOR & POSTERIOR PROBABILITIES 

Prior probabilities are the probabi 1 it ies assigned directly to the 
values of an unknown (random) variable before some particular samp»le 
is taken. Posterior probabi 1 i t ies are the probabi 1 it ies as revised 
in the light of the additional information from that saiTiple- Note: 
this distinction between prior and posterior is always to 

some particular sample. 



RISK 


Usually used inte r changeably with lo&s „ (See also EXPECTED RISK.) 

SAMPLE SPACE 

The space on which the process of s&Q'a&rhttcLl sam.pltri.^ is 
represented- Usually it is an integer space with the x— axis marked 
by the “number of items sampled” and the y—axis marked by “number of 
items found defective”. 

SAMPLING 

The process of randomly selecting and inspecting items from a 
homogeneous lot. 

SEQUENTIAL SAMPLING 

The process of sampling artd inspecting one item at a time until a 
l&r-YTLtri.ctl d.^cis'tc?rL is reached. 

TERMINAL DECISION 

Any of the decisions with which sampltri^ terminates. Usually^ 
"accept lot” and "reject lot" are the only two terminal decisions. 

TYPE I ERROR 

Rejection of the null hypothesis when it is actually true. 

TYPE II ERROR 

Acceptance of the null hypothesis when it is actually false. 



CHAPTER I 


DECISION THEORY AND SEQUENTIAL ACCEPTANCE SAMPLING 

1.0 INTRODUCTION 

The theoretical foundation of (non-Bayes ian ) statistical decision 
theory was laid by Abraham Wald E163, who used minimax procedures 
for minimizing the maximum loss of a decision. Wald also first 
introduced the concept of using decision theory in finding solutions 
to sequential decision problems. We will here be concerned with an 
alternative proposition for making sequential decisions known as the 
B<DLy&£riarL statistical decision theory as explained by Pratt, Raiffa 
and Schlaifer C131. 

In this chapter we introduce the problem of constructing optimal 
sequential sampling plans in the conceptual framework of statistical 
decision theory as originally done by Barnard LZ3 and Wetherill 
E193- Beginning with an introduction to sequential decisions, we 
proceed to a statement of the problem and show how some of these 
plans would appear in practice. We review the salient results 
obtained by earlier researchers working on this problem to bring out 
some of the important issues in developing optimal sequential 
sampling plans. We then present an overview of the present work, as 
developed in Chapters II and III. 
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1.1 UHAT ARE SEQUENTIAL DECISIONS? 

In sequential decision making, statistical decision theory appears 
to be the natural approach for guiding the decision maker toward 
"optimal®* decisions. This may be illustrated as follows. Suppose 
that a decision maker is faced with the problem of selecting a 
decision d from a set of possible decisions D- The effectiveness of 
the decision he makes is dependent upon an unknown parameter p 
(e.g., the quality of the lot being appraised), information about 
which can be ®®pur chased** at a cost of c per *®unit of inf ormat ion" . 
Suppose that the risk in making the decision d when the parameter is 
p is y^(p). The rational decision maker wishes to minimize the total 
risk of his making the decision, which includes the total cost of 
information purchased in making the decision- If he employs the 
Bayesian approach, the decision maker may start with some a. pr-tox't 
knowledge about the parameter p, and at any stage he has the choice 
of either making a ‘'terminal** decision d -e D, or of deferring the 
terminal decision until more information about p is available. Each 
purchase of information about the parameter p represents a **stage*‘» 
The decision maker would proceed sequentially from one stage to the 
lext, evaluate his risk, and terminate the process of decision 
flaking when the total risk of making a decision becomes minimum for 
ill stages. 
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To solve the above problem optimally, i.e,, to develop the sampling 
plan to minimize his total risk it has been suggested that the 
decision maker could use decision theory to guide him in choosing 
the best decision at every stage. A large volume of appropriate 
literature is now available that helps tackle this problem. One may 
refer to DeGroot L31 for a rigorous theoretical presentation of the 
methods. According to Bayesian statistical decision theory, the 
best (i.e- optimal) decision would be the one that minimizes the 
e'Xpe'Ct-eci to toil rtsh t^ie- c:Mrr'&rtt tn/arnuoLtiiorL ahaut pctr'cmu&t&r p. 
The term “risk** here stands for the loss of utility as explained in 
Pratt ^t clI C133. As a special case, risk can represent the 
monetary loss from a decision, provided that monetary loss is a 
satisfactory indicator of the decision maker’s utilities- 

The correct assessment of prior information is as important as the 
effect this information might have on the outcome of the decision 
making process- If a **large** volume of additional information is 
purchased during decision making, the effect of prior information 
will be comparatively low. (How large is *®large** will depend on the 
parameters of the problem.) If, on the other hand, the volume of 
information purchased is *‘low*', the effect of prior information will 
be significant. Barnard £Z1 has argued that the effect of prior 
information will usually be significant in sequential decision 
problems because the objective itself is to minimize the amount of 
information collected, for often the cost of obtaining this 


information is high- 



1.4 


1.2 OPTIMAL SEQUENTIAL SAMPLING PLANS 


Most sequential sampling plans employed in guiding lot 
acceptance/rejection decisions today are those based on Wald’s 
sequential probability ratio test (SPRT) C15I. These plans are 
based on stcitistica.1 measures of consumer’s and producer’s risks 
(expressed as probabilities or odds) at the AQL and LTPD 
respectively. SPRTs suffer from two limitations, however. While 
these plans are very useful in limiting the expected total number of 
items sampled from a lot, they cannot be called "optimal”, because 
they do not minimize any function of cost or risk in the decision 
theoretic sense. Also, they do not offer a generalized model for 
constructing a more accurate representation of the risk to the 
decision maker; SPRTs consider risks specified only at two points — 
the AQL and the LTPD. The above two limitations of existing 
sequential sampling plans can be overcome by constructing plans 
based on a decision theoretic formulation of the problem. 


We first observe that a sequential sampling plan is basically a 
sequential procedure to (optimally) guide decision making about 
accepting or rejecting a lot (or batch) of items that has been 
submitted for acceptance. The two terminal decisions possible here 
are to "accept" or "reject" the lot keeping in perspective the 
unknown parameter, p, known as the lot tyuoLilty. p represents the 



1.5 


true fraction of items that are defective in the lot- Note that p 
is not a measure of the fraction of items that are /oxirL<d defective 
by inspection- In general, these two measures will be different due 
to the presence of inspection errors. The ct priori knowledge of the 
value of p is represented as a prior distribution of p, the unknown 
parameter. Additional information about p may be obtained by 
randomly sampling and inspecting one or more items from the lot, and 
upon inspection classifying them either as defective or 
non—def ective, the cost of each such operation being c. Bayes’' 
theorem is used for updating the prior probabi 1 i t ies based on the 
additional information obtained from sampling- The risk in 
accepting or rejecting the lot when lot quality is p is denoted by 
losses M^Cp) and W^(p> respectively, the objective of the decision 
maker always being minimization of the total expected risk- In 
sequential sampling, after obtaining every sample, the decision 
maker decides optimally on either a terminal decision, or on 
continuing sampling based on his knowledge of the lot quality p at 
that stage of sampling- 

The problem of developing a sequential sampling plan is thus similar 
to the one in Section 1-1 and can be solved optimally. 

To visualize what possibly the sampling plan would be like when the 
above problem is solved, consider a 2~D plane (Figure 1.2.1) with 
the horizontal axis representing the number of items sampled, n, and 
the vertical axis representing the cumulative number of items found 
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defective, r. This description follows the forir^alisms put forth by 
Barnard Cill, CHI and Wetherill C193. The (n, r) space may be 
treated as a representation of the "sample space" of the sampling 
procedure, with each point in it represented as (n, r>- The 
sequential sampling process can be represented as a path CP on the 
sample space beginning at <n = 0, r = 0>, with each successive point 
at stage i having the n~coordinate as i, and the r— coordinate as the 
cumulative number of items found defective after i items have been 
inspected. The entire sample space may be now divided into three 
separate regions, each of which indicates the optimal decision to be 
taken when the path CP enters this region- The three regions 
correspond to the decisions "accept lot," "reject lot," and 
"continue sampling," and thus constitute the sampling plan^ giving 
both the procedure and the criteria for making lot 
acceptance/rejection decisions « 

The boundary between the ctac&pl and eorLtinue regions is called the 
"acceptance" boundary and that between the and c:arLttriXL& 
regions the "rejection" boundary. These boundaries are collectively 
called "decision boundaries". The point at which these three 
decision boundaries meet is called the "meeting point" (Figure 


1 . 2 . 1 ). 
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FIGURE 1.2.1, The sequential sampling process 
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While the procedure for finding 
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the regions for a sequential 
ibed in detail until Chapter III, 
Is of the sequential sampling plan 
to the optimal plan. 


According to Johnson C9II, an optimal sequential 
derived from statistical decision theory comprises 


sampling plan 


the following 

elements s 

(1) Zy The objective of the decision makers 
Usually the minimization of total expected risk. 

(2) Terminal decisions, fd->s 

1 

Usually "accept lot" and "reject lot". 

(3) The Prior distribution ^^(p)5 

Usually that of the lot quality- This distribution can be 
specified as either discrete or continuous depending upon the 
solution methodology. 

C4) Risk, y-(p)5 
D 

The loss of utility (in the decision theoretic sense) associated 
with each terminal decision D at lot quality p. 

i3) The cost of information, cs 

The cost of sampling and inspecting one item. 

(6) The optimal sequential sampling plan, Ss 

This is the resulting optimal sequential sampling plan- 
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These elements of optimal sequential sampling plans and their 
inter-dependencies are shown in Figure 1.E.2. 
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FIGURE Construction of optimal sequential sampling plans- 


A summary of the approaches that variously optimize the sequential 
decision process in acceptance sampling follows. 
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When the prior is represented by a mixed binomial with only two 
componentSy p^ and p^ respectively (which are two possible values of 
the lot quality of a batch of itemsy with r Barnard showed 
that the "most economical" sampling plan (one that minimizes the 
total expected risk) is obtained by a likelihood ratio test 
procedure- If and are the prior probabilities of and p^ 
respectively (a^ a^ = 1 ) , and if ( 1 ~ p^ ) is denoted by q^ y the 
likelihood ratio is defined by 

l(x,y)=X=a^p^q^/a_p^q^ <1.3. 1-1) 

where x defectives and y non-defectives have been observed. The 
risks of errors are given by risks of 
accepting and rejecting the batch at lot quality and p^ 
respectively. One would accept the batch if 


l<x,y) = p^j / pg > a^^/a^ 

and would reject it if 


w % ^ y / X y 

l(x,y) = p^ q^ / p 2 qg 


< a 




<1.3. 1.2) 


(1.3. 1.3) 


while if l(x,y) lies between these two limits, i.e. < X < il- 
would be worthwhile to continue with more inspection. The aim is to 
find X^ and X^ such that the total expected risk is minimized. If we 
denote U by a^^/a^ and L by a^/a^ 2 >^ problem is transformed to 
that of finding U and L optimally. To solve the problem, it is 
required to estimate S(U,L,p), the mean sample size of the 
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sequential procedure, and A(U,L,p) and R<U,L,p) the probabil it ies of 
acceptance and rejection when the lot quality is p. The exact 
formulae from which S, A and R can be calculated have been given by 
Burman C41- It is then possible to derive the equations that 
determine the values for X. ^ and 

1.3-2 VAGHOLKAR AND WETHER ILLS’ 2~P0INT PRIOR HODEL C141. 

Based on basic theory given by Barnard C21, Vagholkar and Wetherill 
have given a method of arriving at the optimum test procedure 
represented by the acceptance and rejection boundaries which for the 
two point prior situation do not meet- CA meeting point is 
desirable in a sequential procedure to rule out the possibility of 
an arbitrarily large sample size before a decision is made.) The 
equations obtained for calculating the decision boundaries are 
similar to those obtained by Barnard (indeed, because the model is 
the same), and are solved numerically by iteration to find the 
values for X^ and X^, The sequential procedure can then be 
implemented using the scoring procedure due to Barnard C13, obtained 
by taking logarithms on both sides of Equation <1-3-1-1)- Two 
illustrative numerical examples are also provided by Vagholkar et - 

As long as the cost of inspection is linearly related to the number 
of items sampled, the above sampling scheme will be optimum- It is 
suggested by these authors that when the cost structure also 
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includes the cost of administrating sequential sampling schemes and 
the cost of pre-processing the items before inspection, the optimal 
scheme may be found by the consideration of group sampling schemes 
having a fixed group size. 

Vagholkar ctl also state that for multi-point prior dist r ibut ions 
with the number of components greater than or equal to three, the 
decision boundaries will in general meet- We feel, however, that 
finding the meeting point analytically may not be easy for k ^ 4. 

1.3.3 WETHERILLS’ CONJUGATE FANILY PRIOR C19, 201. 

In this important work, Wetherill has discussed a number of points 
arising out of the problem of constructing optimum sequential 
sampl ing procedures. 

Wetherill shows that a great simplification in formulating the 
sequential sampling problem is introduced by considering prior 
distributions "closed*® under sampling (i-e., dist r ibut ions whose 
posteriors belong to the same family as the priors with respect to 
sampling from a certain density function — such a family is called 
the "conjugate family"). Suppose that the prior knowledge is 
represented by the distribution of the (vector) decision 
parameter B ^ O, depending on a known (vector) constant a.. Let the 
result of inspection of an item be a random variable X with a 
density function <^(xl6?). After an observation is taken, the 
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posterior distribution of B becomes 

X 

J ^ <^(x. la) da 

o 

If this posterior distribution has the same functional form as that 
□ f ^ (aioi)y say ? (a|/9)y then the prior form ? is said to possess the 
property of being closed under sampling from a distribution - In 
other words, ^ is said to belong to a conjugate family of prior 
distributions « 

^or dist r ibut ions closed under sampling, the sample space of the 
sequential sampling scheme is called the ^—space, the axes of which 
are represented by the vector cx. Thus each member of the family of 
dist r ibut ions ^(a|ot) is represented by a point on the space. 
Hetherill derives the equations which lead to the meeting point on 
the decision boundaries on the ^-space, and gives a dynamic 
programming scheme to find the decision boundaries starting from 
this meeting point- 

The equations leading to the meeting point are obtained as 

f cy iB} - y (Byj da = o (1.3.3.1> 

r a 

and 


f qi€A,&) CW W) - W de 

r a 


o 


( 1 . 3 . 3 . 2 ) 
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where, W (d) and W (6?) are the risks of making decisions accept and 
reject respectively given that the decision parameter is q (ca ) 
is the probability that for an a on the meeting point boundary^ one 
observation leads to points in the acceptance region. 


Wetherill proves that both the decision boundaries do in fact meet 
at the meeting point obtained from Equations <1.3. 3.1) and 
(1.3. 3-3). The optimal t&rminjal risk at any point ft in the space 
is given by 


R^(^) 


= min 




(1.3. 3. 3) 


while the cont ion. risk can be calculated optimally by 


R^(/9> = 1 + JJ ?(6'tri> <j!.(xl<9) dx dt» 


(1.3. 3. 4) 


where ft corresponds to the posterior distribution for a prior 
distribution ft and an observation x. The optimal risk at ft is then 

Rifty = min ^ R^(/9>, (1.3. 3. 5) 

Now, if point ft is reached after n sampling steps, the risk R^(/5) 
cannot be calculated unless the optimal risks of all points ft^ 
reached in (n + 1) steps is known. If N is the sample size 
corresponding to the meeting point, then decisions at all points 
reached in N steps (one such point is the meeting point) have to be 
terminal decisions, and the minimum risk at each such point can be 
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calculated easily as in Equation Cl. 3. 3. 3). Thus it is possible to 
locate the meeting point at CM, R) by solving Equations Cl. 3-3.1) 
and (1-3-3. 2) simultaneously, and then work backwards in n to find 
the optimal risks of all points co r r es ponding to sample size (N—l), 
(N-2), n, O, using the Equations (1.3.3. 4) and (1-3-3-5)- 

The above principle gives an optimal sampling plan, because it 
conforms to the principle of optimality given by Bellman C33- 

1-4 THE PRESENT WORKS THE K-POINT PRIOR WITH K > 2. 

In this work, we find optimal sequential sampling plans based on the 
general k— point prior distribution with k > 2 using the concept of 
the meeting point developed by Wetherill E191, and thus extend 
(1) the range of priors available for developing sequential plans by 
including the k-point prior distributions, and, (2) the ability to 
work with prior distr ibut ions with the number of parameters greater 
than 2- 

We study the following issues in the construction of optimal 
sequential sampling plans using the k-point form of the prior 
dist r ibut ion s 

1) Formulation of the problem of obtaining optimal sequential 
sampling plans and the derivation of equations for 
obtaining the decision boundaries. 
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E) Computational aspects of the construction of the decision 
boundaries for the optimal sequential sampling plans. 

3) Calculation of the ASN and the OC curves for the k~point 
prior case- 

4) Calculation of the expected risk in using the sampling 
plan. 

5) Sensitivity of the plan parameters Ce-g., risk, ASN, Pa> 
to the following factors^ 

1) The cost of inspection. 

2) The parameters of the prior distribution 

3) Location of the meeting point. 

4) Imperfect knowledge of the inspection cost 

3) Imperfect knowledge of the decision losses 

6) Imperfect inspection. 

6) Exact analytical solution for the meeting point for the 

2~ and 3-~point prior case. No meeting point is found for 
the k = 2 case but a solution is obtained fay forcing a 
meeting point at a sufficiently large n. This is then an 
alternate way to approach the Vagholkar clI C143 


solution based on the likelihood ratio test. 



CHAPTER II 


FORMULATION OF OPTIMAL SEQUENTIAL SAMPLING PLANS BASED ON 
THE K-POINT PRIOR DISTRIBUTION 

2.0 INTRODUCTION 


The most popular cantinuoxiS' prior dist r ibut ions used in the 
literature for deriving optimal sequential sampling schemes have 
been the beta and the gamma dist ribut ions , both specifiable 
completely by two parameters. Lindley and Barnett E12I1 have given 
sequential plans and also numerical results for the beta prior 
distribution, sampling being done from a binomial population, 
assuming the losses of the terminal decisions to be linear. Plans 
for the gamma-Foisson case have been given by Lechner E113- The 
diiscr-ete prior distribution for which analytical results are 
available is the 2— point prior dist r ibut ion . Vagholkar a.1 C143 
derive the equations for the 2~point prior case using the likelihood 
ratio test given by Barnard C23, and outline an iterative procedure 
for obtaining the decision boundaries- 


In the present work, we derive optimal sequential sampling schemes 
using the generalized k— point prior distribution assuming arbitrary 
loss functions. We use the meeting point approach suggested by 
Wetherill E191 to find the decision boundaries. 
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As pointed out in the literature, a k— point prior distribution 
offers the following advantages over the more commonly used 
conjugate prior distributions? 


1 ) 

The k— point 

prior 

is a less 

restr i ct i ve 

form 

for 


representing 

prior 

knowledge . 

Other conjugate 

prior 


distr ibut ions 

( such 

as beta, gamma 

and normal ) 

have 

only 


two parameters and have been argued against by Barnard C2II 
who showed that prior dist r ibut ions based on only two 
parameters might give inconsistent results and commented 
that such dist ribut ions are not a representation 

of prior knowledge. The k—point distribution has (Ek - 1) 
parameters that may be independent ly adjusted to best 
represent prior information- 

E) Barnard CEIl also points out that in real life, the 
distribution of fraction defective in a lot will usually 
be a multi— point distribution, especially when the lot 
comes from a process operating at several different levels 
of quality- 

For k = 3, the resulting model can be solved analytically and 
exactly to produce the meeting point. As we show, for k 4, the 

model can be solved numerically to converge to the exact location of 
the meeting point of the acceptance and rejection boundaries. ye 
find, however, that the numerical procedure does not always converge 
to a solution (if it exists). This may require further study. If 
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the meeting point is forced at a sufficiently large H, the optimum 
decision boundaries can be found optimally even if the meeting point 
cannot be found analytically or numerically. 


2.1 FORMULATION OF THE k-POINT PRIOR MODEL 


In the k-point prior model:, it is assumed that the unknown lot 
quality (fraction defective) p can have k distinct values p^, p^y 
p^ - yith each possible p^ , the decision maker associates a 
prior probability of its occurrence, a^ , based on the knowledge 
existing about the possible distribution of lot quality. Thus the 
k“point prior distribution P of the lot quality p becomes 


PCp = p^ ) = a^ for i = 1, 2, k. 


i=k 

where ^ a^ = 1, k 2- (Figure 2.1.1) 
i = 1 


^k-1 


Probability j 


( . 


Pe P3 


pk-1 Pk 


lot quality (p) *■ 


(E.1. 1) 


FIGURE E.1.1. The k-point prior distribution. 
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Sequential sampling to guide acceptance/rejection of the lot would 
proceed by sampling and inspecting one item at a time at every stage 
of sampling, the cost of this operation being c each time an item is 
inspected (c is assumed to be known to the decision maker). The 
decision at each stage would be guided by a sampling plan, which 
would be represented by the acceptance, rejection and continuation 
regions on the sample space as discussed in Section 1.2. Let the 
three decisions possible at every stage, namely "accept lot", 
"reject lot", and "continue sampling", be abbreviated by notations 
"ac", "rj" and "cn" respectively. The losses associated with each 
of these decisions for any given lot quality p may be denoted by 
W^j.(p), W^j(p) and W^^(p) respectively and are assumed to be known 
to the decision maker. The objective of the decision maker is to 
sample items one by one to eventually accept ar reject the lot such 
that the total (expected, Bayes’) risk is minimized. 

We will attempt to find the optimal decision regions for the above 
sequential decision problem following the procedure given by 
Wetherill C193. The procedure begins by locating the meeting point 
of the decision boundaries, and subsequently uses an iterative 
scheme (a recursive scheme converted into an iterative scheme) to 
find the three decision regions for all sample sizes less than or 
equal to the maximum possible sample size (the maximum sample size 
corresponds to the sample size required to reach the meeting point). 
This procedure results in optimal decisions (decisions which 
minimize the total expected risk) if the meeting point can be 
accurately located (Wetherill C193). 
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2.H EQUATIONS FOR DETERMINING THE MEETING POINT 


At the meeting point (n , r ) of the decision boundaries, by 

mm 

definitioHr the risk associated with sampling one more item is equal 
to the risk of making either of the two terminal decisions^ namely, 
accept or reject the lot. Also, the risk of continuing with 
sampling is equal to the cost of taking the next sample plus the 
risk of making an optimal terminal decision subsequently. If q^^. (p) 
denotes the probability of rejecting the lot at the meeting point, 
the risk of continuation at the meeting point is 

R (p) = q .(p) y ,(p) 4- (1 - q .(p)) y Cp) C2.2.1) 

c ^ ^rj ^ rj ^ ^rj ^ ac ^ 


The tata.1 €?y:p^ct&d. rt&h of continuation at the meeting point is 

found by “weighting" R^(p) with the k-point posterior distribution 

at (n , r ), and adding the cost of sampling c, 
m m 


i=k 


R = c 
cn 


y q .Cp.) y .(p.) 4- ( 1 - q .Cp.)> y Cp.) 
L j_ rj rj ^rj ac 


a - (n , r ) 
1 m m 


i=1 


<2.E.E) 


r n -r 

, vmm 

a. (p.) C1— p.) 

1 1 1 


where ^,*Cn_,r^) = 
1 m m 


i=k 

2 a, (pj)'""' (1 - pj) 

i = 1 


n — r 
m m 


(2.2.3) 


.vp J J 
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is the posterior distribution at (n^ ^ ' obtained by applying 
Bayes Theorem, having started the sampling process with the prior 
(2.1.1). (cf page 282, yetherill C201.) 


The expected risks of the terminal decisions accept and reject at 
the meeting point are given by 


ac 


i=k 

= 2 

i = 1 


i=k 


W (p. ) 
ac ^1 


a - (n , r ) 
1 m m 


and R 


rj 


U -(p-> a.Cn ,r ) 
rj 1 mm 


(2.2.4) 


i=1 


At the meeting point, R = R . = R , and q .(p) = p. Equating R 
^ ^ ac rj cn ^rj ^ 

and R , and substituting for q (p), 
cn ^ ^rj ^ 


i=k 


c + 


2 


w . ( p. > + 

rj 


(1 - 


p. >y (p 

ac 


i=1 




a. (n , r ) 
1 m m 


i=k 

= 2 

i=1 


W (p.) a.(n ,r ) 
a c 1 1 m m 


or. 


i=k 


y y (p.) “ y .(p.) - c/p. a.(n ,r ) = 
Z [ ac rj J 1 m m 


i = 1 


(2.2.5) 


Because we are at the meeting point, we can equate R and R 

^ ^ ^ ac rj 

resulting in 


2 [ 

= -! 


y 

ac 


CPi> 


- W (p. ) 
r 1 


a . ( n 5 
1 m 


r ) = 
m 


O 


( 2 . 2 . 6 ) 
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Substituting for in equations (2.2.5) and (2.2.6), we get 


i=k 

2 [ 

i = 1 


W (p. ) 
ac 


W -(p.) 
rj 


c/p. 


p.a.Cp.)*^* (1 — P-) 

1 


n -r 
m m 


= O (2.2.7) 


2 [ “ac'Pl’ - «rJ<Pi>] 'fi''" ^ 

i = 1 


Pi> 


n — r 

IT) II) 


0 


( 2 . 2 . 8 ) 


Equations (2.2.7) and (2.2.8) are two simultaneous non-linear 
equations in two variables (n , r^) — the coordinates of the 
meeting point. If the meeting point exists (conditions not yet made 
clear in the literature), solving these equations will in general 
be possible numerically using any of the various quardrature methods 
available for solving such equations in two variables. The 
Newton— Raphson procedure is one possible approach which could be 
used here: we have found that it is relatively simple to implement 
and converges quickly. For the description of this procedure, one 
is referred to Kelly C10I1, or to any other standard text on 
numerical methods. Analytical solutions are possible for small 
values of k — the exact solutions for k = 2 and k = 3 are presented 


in Chapter II. 
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2.2.1 INTERPRETATION OF THE MEETING POINT 


It is necessary here to clarify the concept of the “meeting point” 
in the context of the sample space (Section 1.2). ye will use here 
the following standard notations: 


c&tlix} : the smallest integer greater than or equal to x, and, 
/loorCx): the largest integer less than or equal to x, 
where x is a real number. 

The sample space is a discrete space, with both n and r taking only 

integer values. The solution of Equations (2.2.7) and (2.2.8), if 

it exists, in general will yield non-integer (real) values for n^ 

and r . Hence the “real” meeting point (n , r^) is not realizable 
m ^ ^ m m 

during sampling. But these (real) values of n and r give us two 
^ ^ ^ mm 

trivial, but important results: 

1) Acceptance is the optimum decision at all points 

(C4etl(n ),r ) and rejection is the optimum decision at 
m a 

all points (e<eil(n ), r ), where r < /loorCr ) and r > 
^ m r a m r 

c&tlir ) (see Figure 2. 2.1.1). 
m 

2) The maximum possible sample size N is ee-il(n^), beyond 
which only terminal decisions (i.e, “accept” and “reject”, 
but not "continue” sampling) are optimal- 



(N, c«.iL(r )): "REJECTION" OPT; 

“ m 

(N, /ioor(r )): "ACCEPTANCE" OF 

FIGURE 2.2.1. 1. Interpretat ion of the meeting point- 

The above two results will be used in finding the optimal d 
regions which minimize Bayes* risk in the following s 

Henceforth^ we will denote e<eil(n ), the sample size at the 

m ^ 

point, by N, the maximum sample size obtainable. 

2.3 CONSTRUCTION OF THE DECISION REGIONS 

In this section we derive an iterative procedure fo 

construction of the decision regions for an optimal seqi 

sampling plan. At any point (n, r) on the sample space, the 

of the three possible decisions are given by R (n,r), R .(n,) 
^ ^ ac rj 

R (n,r) where subscripts ac, rj and cn have their usual me^ 
cn 

The optimum risk at this point R (n,r> is given by 

R*(n,r) = min - R (n,r), R .<n,r>, R (n,r) v 

a c • J cn 


V 



I 


?re R (n,r) and R .(n,r) can be calculated as 

« c r* j 

i=k 

_(n,r) =) W (p.) a.(n,r) and 
c Z ac 1 ' 

i = 1 

i=k 

.(n,r) =yw .(p.) a.(n,r) <2.: 

j Z rj ^1 1 

i=1 

d R (n,r) is obtained by the recursive relation 
cn ^ 

(n,r) = c + -I q .(n,r) R*(n+1, r+1) + (1 — q .(n,r)) R*(n+1, r 
n rj ' rj 

( 2 .: 

ere q^j(n,r) is the probability that the next (n+1-th) 
mpled will be found defective and hence the lot rejected at 
age of sampling. 


(n+1, r+1) 



FIGURE 2.3.1. The result of inspecting one item at (n, r) 

defective item will lead to the point (n+Irr+D on the sai 
■ace while a non-defective item will lead to the point (n+1 
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Thus if the decision taken after taking one more sample at (n,r) is 
an optimal decision^ the risk R^^(n,r) is obtained as in Equation 
(2-3-3). In this Equation (2-3. 3), q^j(n,r) is calculated from the 
posterior distribution as 

i=k 

qi^j(n,r)=2p^a^(T»,r) <E.3.4) 

i = 1 

The optimal decisions for all r at N are terminal decisions, the 
risks of which can be calculated from the Equation (2.3.2). Thus 


ac 

(N, 

r ) 

for 

r < /laoT' ( r ) 
m 


(N, 

r ) 

for 

r > 1 i r ) 

m 


(2.3.5) 


A recursive scheme based on Equation (2.3.3) will terminate at the 
N-th stage, where R (N,r) will be obtained from equation (2.3.5) 
without further recursion. Every (n, r) point in the space will 
thus become marked as an acceptance, reiection, or a continuation 
point, depending upon which risk is found to be minimum. 
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2.3.1 SOME COMPUTATIONAL ASPECTS OF THE CONSTRUCTION OF DECISION 
REGIONS 

2.3.1a Conversion of the recursive scheme to an iterative scheme. 

The recursion scheme (2.3.3) starting at (n = O, r = O) will require 
a large amount of storage (in the form of a recursion stack) and 
will usually be slower. Because it is known beforehand that this 
scheme will terminate after exactly N steps, it can be converted 
into an iterative scheme. This iterative scheme begins with 
calculation of the optimal risks at sample size N from Equation 
(2.3.5) and proceeds bctckxMir-d^ in n to find the optimum risks at all 
the sample points with 0 n ^ N. 

2.3.1b Computational effort. 

The most expensive computational element in the above scheme is the 
calculation of the posterior distribution at every point (n,r). 
Enumeration of risks at all points with n N in the sample space 
will require that the posterior distribution be calculated at 
(N+1 ) (N+E ) /2— 1 points (as points with r > n cannot be reached). But 
it is possible to restrict this enumeration to far fewer points in 
practice because points "beyond" the decision boundaries will never 
be reached during the sampling process (sampling will terminate with 
a decision as soon as a boundary is reached). If w(n) denotes the 
number of points in the continuation region for any n, and the 
average "width" of the continuation region ( w(n) averaged over all 
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n < N) is w, the number of points at which the risks have to be 
enumerated are (w + 2)N, which is the sum of all continuation and 
decision boundary points. 

(We have found w to be usually a small number less than 10 in most 
problems. In such cases, the posterior needs to be calculated at 
0(N) points instead of 0(N ) points. w is largely a function of the 
magnitude of the risks, and is large when the risks are large.) 

2.3.1c THE LARGEST POSSIBLE "EFFECTIVE” SAMPLE SIZE, N* . 

Sequential sampling will always terminate at or before a certain 

•Jt* 

sample size N , defined as the minimum sample size for which there 

^ -If- 

is no continuation point at (N -^1). In general, N will be less 

than or equal to N (N = c&ilin )), and we illustrate this using an 

m 

hypothetical example of decision boundaries shown in Figure 2.3. 1.1. 
Assume that the decision boundaries can be approximated to 
continuous curves (shown as straight lines in the figure) which 
separate the acceptance, the rejection and the continuation region. 
In the illustration, in the vicinity of the (exact) meeting point, 
the decision boundaries are relatively "narrow" and do not contain a 
meeting point between them at sample sizes (N-1>, (N-3) , (N— 4) and 

(N-7). However, there are continuation points at all sample sizes 
less than (N-7). Thus, (N-7) is the maximum obtainable sample size, 
N*, for the sequential sampling plan represented by boundaries shown 
in the figure. There will, in general, always be a continuation 
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E ■ ■ 

j 

I . 

N-7 N-5 N-3 N-1 N 

Number of items sampled (n) 


FIGURE 2. 3. 1.1. The effective sample size, N*. 


point between the two decision boundaries if the separation between 
them is greater than one unit on the vertical scale. 


2.4 "BOUNDS" : THE ALGORITHM TO COMPUTE THE DECISION BOUNDARIES. 

Algorithm BOUNDS, given below may be used to compute the decision 
boundaries for a k-point prior model. This algorithm is implemented 
in C as shown in Appendix A. Taking the analytical meeting point, 
the losses, and the prior distribution as the input, BOUNDS 
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calculates the effective maximum sample size, N*, and the 

acceptance and rejection boundaries on the decision space. 

The inputs to BOUNDS are: 

1> The meeting point, (N, R), where N=eeil(n ) and 

m 

R=/loor(r >. Recall that (n , r ) is the analytical meeting 
m m m ^ ^ 

point . 

E) The risks of accepting the lot W (p) and rejecting the lot 
W^j{p) as functions of the lot quality p. 

3) The prior distribution of p (Equation 2.1.1). 

4) The cost of sampling per item, c. 


The output of BOUNDS are: 

1) The effective sample size, N . 

E) Acceptance and rejection boundaries for all sample sizes n 
i N , which are obtained as arrays D^Cn3 and D^^CnJ 


respectively. 
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algorithm BOUNDS 

BEGIN 

* 

N ■« — N; /* initialize effective sample size to 

MAX_R 4 — inin(N, R + 1); /* upper bound on the cumulative number 

defectives for sample sizes less than ! 
D^CNO 4 — MAX_R-1; /* acceptance boundary at N */ 

DpCNU 4 — MAX_R; /* rejection boundary at N */ 

tor r 4 — 0 -to D^CND do 

i=k 

R CN, r3 4 — ^ ^ac^^i^ posterior(N, r, i); 

i=1 

end 

/* above calculates the acceptance risks (optimum risks) at 
acceptance points at N */ 

tor r 4 — R + 1 to R + 2 do 

i=k 

R CN, R+m 4 — T W .(p.) posterior(N, R+1, i); 

Zj f" J 1 

i = 1 

end 

/* calculate the risk of rejection (optimum risk) at (N, 
and (N, R+E) */ 


N */ 
of 
I */ 


all 


R+1) 


/* backward calculation of risks begins henceforth */ 
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for n 4 — N - 1 to O do /* n indeKes the sample size */ 


for r 4 — O to MAX_R do /•» r indexes the number of defectives */ 


y W (p.) posterior prob (n, r, i); 

Zi a I- 1 


/* acceptance risk at (n,r) */ 


R_- 4 — ) W .(p ) poster ior_prob (n, r, i) 
i = 1 

/* rejection risk at (n,r) •»/ 


^ 2 [*^i RCs+1DCtDj X 

i=1 poster ior_prob (n, r, i); 

/* the continuation risk at (n,r). c is the cost of sampling 
andq. =(1— p.) */ 


R Cn, rl 4— min •{ R , R R 

ac rj cn J 

/* the minimum risk at (n, r) */ 


If R Cn, r3 = R then d 

ac ac 


d + 1 ; 
ac 


/* d^^ counts acceptance points at n ■»/ 


if R Cn, r3 = R then d 

cn cn 


d + 1 ; 
cn 


/* d counts continuation points at n */ 
cn 


end 


/# the loop for r */ 
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.CnI 4 — d ; 
4 ac 


,CnI <— d 
? ac 

* calculate 


+ d 
the 


+ 1 ; 
cn 

acceptance 


and 


rejection 


boundaries at n */ 


^ — DpCnD; /* new upper bound on the cumulative number of 
defectives for sample sizes less than n */ 

f d =0 then N* 4 — n ; 
cn 

/* if there are no continuation points at n, 
it becomes the new meeting point. See Section 2.3.1c ■»/ 

nd /■» of the n loop */ 


ND. /* of algorithm BOUNDARIES */ 


!.5 THE ASN AND THE OC CURVES 

Algorithm "BOUNDS" calculates the decision regions for an optimal 
sequential sampling plan. We will now discuss the procedure for 
:alculating the average sampling number (ASN) and the operating 
characteristic (OC) curves for the present sampling plan. The ASN 
is defined as the expected number of items that are sampled before a 
terminal decision is reached. The OC curve gives the probability of 
accepting the lot at any given lot quality p. As we shall see, it 
is possible to numerically compute the ASN and the OC curves using 
our procedure foecct-use the continuation region for the sampling plan 
is bounded as sampling would not need to beyond the maximum 
effective sample size N . 
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E.5.1 "REACHABLE" AND "NON-REACHABLE" POINTS 

Points on. the decision hoxindax-ies belong to two mutually exclusive 
(and exhaustive) sets: one consists of "reachable points" and the 
other of "non-reachable points". All reachable points on the 
decision boundaries are necessarily termination points. 
Non-reachable points are so called because they cannot be reached 
during any sampling process, being "blocked from reach" by other 
decision boundary points. 

We will illustrate the concept of "reachability" and 

"non— reachability" of points with the help of an example. Let (n, 

r^^(n)) be a point on the acceptance boundary and (n, r^j(n)) be a 

point on the rejection boundary at some stage of sampling n. In 

Figure S.5.1, the acceptance point r (6) is a non-reachable point 

because it can only be reached from points (5,0) and (5,1), both of 

which are points at which acceptance will be preferred, and sampling 

will stop (in fact, point (5,0) will never be reached). The 

acceptance point *"^^(7) is non-reachable because it can only be 

reached from (6,0), which lies in the acceptance region, and (6,1), 

a non-reachable point. Similar arguments can be given for points 

r .(9) and r .(10). 
rj rj 

The criteria for a point to be non-reachable is: 

1) r (n) is a non-reachable point iff r (n) = r (n— 1). 
ac ac ac 

r .(n) is a non-reachable point iff r_^.(n> > r .(n-1). 
rj rj rj 


2 ) 
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0 1 2 3 4 5 i 7 « 10 11 12 13 14 15 

number of Items sampled (n) 


FIGURE 2.5.1. Reachable and non-reachable points. 

2.5.2 PROCEDURE FOR THE CALCULATION OF OC AND ASN CURVES 

A sampling process can terminate only at a decision boundary. 

P (n,r) denotes the probability of reaching any point (n,r) on 
R 

sample space with lot guality p (Eguation 2. 5. 2. 5), then 
probabilities of termination of sampling at a point on the deci 
boundary for sample size n are ^R^"' r^^Cn)), 


If 

the 

the 

sion 


for 
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a given lot quality p. Then the probability of acceptance of the 
lot, i.e., the OC curve is given by 

n=N 

Pa(p> = J ""ac^"^^ (2. 5. 2.1) 

n=1 

and the average number of items sampled at lot quality p is 

n=N 

ASN{p) = J ] (2. 5. 2. 2) 

n=1 


where it should be noted that if r (n) or r .(n) are non— reachable , 

ac r j 

we have: 


P„(n, r (n)) =0. 
R ac ^ 

and j 

P|^(n, r^j. (n) ) = 0 " 


( 2 . 5 . 2 . 3 ) 


We now turn to the determination of general. The 
probability of reaching cLny point (n,r) in the sample space by 
starting sampling at (0,0) is equal to the sum of the probabilities 
(each sampling being independent) of following each of the unique 
paths that lead to and end at that point. The probability of 
reaching any point in the terminal decision region which is not a 
boundary point is O. Only the points in the continuation region and 
the points on the decision boundaries are potentially reachable. 
Thus, we can calculate Pp(n,r) for all points (n < N, r < n) in the 
decision region using the following recurrence relation: 
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Pp(n,r> = p Pj^(n-1, r-1) + (1 - p) Pj^(n-1, r) (2. 5. 2. 4) 

with the following boundary conditions: 

Pj^(0,0) = 1 (2.5.2.5a) 

PD<nrry = O for r St r .(n) + l and r :£ r (n)-1, and, (2.5.2.5b) 
^ r j 3 c 

P„(n, r (n))=0- 
R ac ^ 

and • for non-reachable boundary points (2.5.2.5c) 

Pp, (n,r.(n))=0-' 

R rj 

Sequential sampling will eventually always terminate at a point that 
lies on one of the decision boundaries. One way of looking at this 
is the following. Consider the following transition matrix, 

constructed using the Equation (2. 5. 2. 4) and Conditions 

(2. 5. 2. 5 a, b, and c), showing the tion proba.biliti&£ on the 

sample space when ort& item is inspected 


to: 

(0,0) 

(0,1) 

. . 

(n y n ) - - 

( n+1 , r ) 

(n+1. 

r+1) .. 

* 

(N +1 

.y) 


/ r-om: 






- 





(n , r ) ’ 

0 

0 

.0. 

1 

o 

• 

o 

1-p 

P 

.0. 

0 

* ROW 

1 

(n,r (n)) 
ac 

0 

0 

.0. 

■ 

O 

• 

O 

0 

0 

.0. 

0 

ROW 

2 

(n , r . (n ) > 
rj 

0 

0 

.0. 

o 

b 

0 

0 

.0. 

0 

ROW 

3 

( N* , X ) 

0 

0 

.0. 

0 .0. 

0 

0 

.0. 

0 

ROW 

4 


FIGURE 2.5.2. 1. Transition probabilities on the sample space. 


i 
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The transition probability from any boundary point (ROWS 2 & 3) is 
always zero, and sampling will thus stop at a boundary point. Also, 
no transition is possible once the effective sample size N is 

reached (ROW 4>. Because a “reverse" transition (i-e-, (n,>c) ► 

(n-1,y)) occurs with a zero probability, sequential sampling 
beginning at (0,0) will result in a (strictly) monoton i cal ly 
increasing sample size and will terminate irt th& worst coso at the 

■jf- 

effective sample size N (from which no transition is possible), if 
a decision boundary point is never encountered during the course of 
sampl ing . 


2-6 EXPECTED RISK OF THE OPTIMIZED SEQUENTIAL SAMPLING PLAN 


The “ideal** decision (based on perfect knowledge of p) will always 
have a risk that equals the minimum of W (p) and W .(p) (see 
Section 2-1) at all values of the lot quality p- The optimum risk of 
such a decision is denoted by 


R*(p) = 


min ^ 


W (p) 
ac ^ 


W -(p) 
rj 


} 


( 2 . 6 . 1 ) 


■a- 

Where acceptance will be preferred when R (p) = W (p) and rejected 

o C- 

when R*(p) = W (p). 

ac 


We assume that it is not economically feasible to make the ideal 
decision because p is an unknown and the cost of obtaining 
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information about p may be high. The "economically optimum" 
sampling plans obtained using the results of the previous sections 
will in general have an expected risk R(p) > R*(p). The expression 
for R(p) can be obtained using Pa(p) as follows! 

R(p) = Pa(p) W (p) + (1 - Pa) W .(p) (2.6.2) 

a.c r J 

with Pa(p) being calculated from Equation (2. 5. 2.1). 

The shift of R(p) from the ideal is the "loss due to imperfect 
information" and is defined as: 

& = R(p) - R*(p) (2.6.3) 

It is easily seen that an optimization of R(p) subject to the cost 
constraints of sampling is equivalent to a minimization of 6, 
because R (p) is a constant dependent only on the actual lot quality 
p. The optimum sequential sampling plan also minimizes <5 . 

In general, there can be many factors that could affect the 
optimality of the sampling plan. Some of these factors are 
imperfect inspection, imprecise knowledge of the various risk and 
cost parameters, and the subjective differences in the assessment of 
the prior distribution. The location of the meeting point is another 
factor which might affect the plan’s optimality — when the meeting 
point cannot be exactly determined or does not exist and may be 
forced arbitrarily at some large enough sample size to stop the 
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sampling process. Some of these factors might even lead to a 
in 6 as we shall see during sensitivity analysis! The 
following section is devoted entirely to one such factor, the effect 
of imperfect inspection (as items are sampled one by one) on The 
optimality of the plan. 


2.7 THE EFFECT OF IMPERFECT INSPECTION 

Two types of errors are possible during inspections rejection of a 
non— defective item (Type I error) and acceptance of a defective item 
(Type II error). Studies show that the above two types of error can 
at times adversely affect the optimality of decision making during 
acceptance sampling (see Dorris a.1 ) . 

If the probability of the Type I error in inspection is and that 
of the Type II error is €>, the fraction of the items that will be 
classified as "defective" (known as the apparent fraction defective) 
will be 

p = p (1 - e) + (1 - (2.7.1) 

If these & and ip errors are not corrected, the decisions regarding 
the acceptance and rejection of the lot will be based on this 
apparent lot quality p instead of on the true lot quality p. the 
extra loss incurred due to imperfect inspection will then be 
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<5 - R(p) - R(p) (2.7.2) 

If the cost of correcting imperfect inspection from levels ©) to 

levels (^, &) is given by c(<^, &} , the gain resulting from 

this correction is given by 

G = R<'p) - R(p( 1 - 6 ») + ( 1 -p)<^) - clip, e, ip, &) (2.7.3) 

When p is known, or can be estimated, a correction lip, — *■ lip, 6) 

should be applied if G > 0. 


p 

PROCESS » INSPECTION 

fraction • 
defective 


p(1-6i)+<t(1-p) 

fraction cctlLi^d. 
def active 


pd + ( 1-p) ( 1-1^ ) 

fraction call&d. non— def active 


FIGURE 2.7.1. The effect of inspection error. 


It may be possible to derive the optimal sequential sampling plan 
given IB ip). However, we will not attempt that here. 



CHAPTER III 


ANALYTICAL SOLUTIONS FOR SOME K-POINT PRIOR MODELS 


3.0 INTRODUCTION 

While the k— point prior case will in general be usually solved 
rtxmu&r-icctlly to yield the meeting point and the subsequent 
determination of the decision boundaries, two particular cases, 
those of k = 2 and k = 3, can be investigated . The 
discussion of the 2-point prior case can also be found in Wetherill 
C20D and Vagholkar ctl C14I. Wetherill shows that the optimal plan 
for the 2— point prior case is an SPRT, and that the meeting point 
does not exist when the unit cost of sampling is a constant. 
Vagholkar et a.1 use the likelihood ratio test procedure to find the 
decision boundaries for the 2— point prior case. 

Here we discuss the 2-point prior case using our formulation (which 
is similar to that of Wetherill C193) and also present an exact 
analytical solution for the meeting point for the 3-point prior 


case . 
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3-1 SOLUTION FOR THE 2-POINT PRIOR DISTRIBUTION CASE 


The 2— point prior case is the simplest of all possible k-point 
distributions. The analytical solution for this case based on the 
likelihood ratio test of Barnard C2D was give by Vagholkar &t clI 
C14I1. (For the 2— point prior case, it was shown by Barnard that a 
likelihood ratio test can be constructed to produce the expected 
minimum cost sampling plans, which he called optimum plans.) 


Using the Equations (2.2.7) and (2.2.8) we show in this section that 
the meeting point does not exist for this case. These Equations 
reduce to the following form: 


i=2 


y ' W (p.) - W .(p.) — c/p. p. a. (p.) "^ ( 1 - p.) 
L ac rj 1 *^1 1 


n — r 
m m 


i = 1 


= 0 

(3.1.1) 


i=2 


Y r Ul (p.) - W (p.)] a. (p.) (1 - p.) 

Z L ac rj J 1 


n -r 
m m 


(3.1.2) 


i=1 


On substituting 


1 ^ Kc'Pe' - “rj'Pg’ - 


-1 


(1 - p^) and 
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<1 - P2> 


in (3.1.1) and (3.1.2) and taking logarithms we obtain two 

simultaneous linear equations in n and r , as follows. 

m m 


r In 
m 

[p^qg / Pgq, 

+ In q, / qg] 

r In 
m 

P^qg / Pgq,] 

"m *"[‘’1 ' Xs] 


ln[agkg / aj 


(3.1.3) 


(3.1.4) 


The above two equations have no solution, because the determinant is 
0. Therefore a meeting point cannot be found for the 2— point prior 
distribution. If the cost of sampling c is dependent on n^, 
then k^will involve n^, and the Equations (3.1.3) and (3.1.4) 
will in general have a solution. 

Wetherill C203 points out that a meeting point for the 2— point prior 
distribution will in general exist if the cost of sampling, c, is a 
function of n — the number of items sampled. Wetherill finds the 
meeting point when c(n) = n , and shows that if the meeting point 
occurs at sample size N, N goes to infinity as y — ^ O- 

In many practical situations, c will be a constant and therefore a 
meeting point at a finite N will not be found. In such situations, 
the decision boundaries may be obtained by the likelihood test ratio 
procedure given by Vagholkar clI C14I1 will have to be used. 
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Wetherill CH03 also suggests that an approximate solution for the 
decision boundaries may be found by forcing the meeting point to 
occur at a "large enough" sample size. We feel that this will 
require that a heuristic be developed for estimating what the value 
for this "large" sample size n^ should be. We suggest below one 
such heuristic for finding a "good" meeting point '"m^' 

Equation (3.1.3) gives a straight line (in the decision space) on 

which the decision maker is indifferent between acceptance or 

rejection. Thus to be a meeting point, (n , r ) must lie on this 

mm 

line. Hence for any given n there must exist a unique r . If 

R(p,n^) represents the risk of a sampling plan obtained using 

(n r ) as a meeting point when the lot quality is p, then the 

following (generic) algorithm will produce a close enough estimate 

f or n ! 
m 



III. 5 


heuristic FIND MEET PT 


/•* -finds "meeting point" sample size n such that the risk of the 

m 

sampling plan is relatively insensitive to n */ 

m 


BEGIN 


n^— O 

repea-t n — n + STEP 


/■» STEP is a user-specified positive integer ■»/ 
(p, n) - R(p, 


un-tl 1 max 
P 


(0!<pt<1> 


/* where wc is a small positive real number */ 



n . 


END. /*ofalgorithmFIND_MEET_PT«/ 


Once a meeting point is forced at (n , r ) , numerical procedures as 

mm 

given in Sections (2.3), (2.4), (2.5) and (2.6) can be used to 

develop the decision boundaries, and hence the sequential sampling 
plan. We will derive a better technique for forcing the meeting 
point at a sufficiently large (n ,r ) in Section 4.4 (3), when we 

discuss the sensitivity of the sampling plan to the location o-f the 


meeting point. 
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3.2 EXACT SOLUTION FOR THE 3-POINT PRIOR CASE 

The prior distribution in the 3— point prior model is specified as 

i=3 

for i = 1, 2 and 3, where ^ ^i ~ 

i = 1 

Here , Pg and p^ are the three possible values of lot quality. 


probabi 1 ity 


Pl P3 


FIGURE 3.2.1. The 3— point prior distribution. 


The 3-~point prior distribution appears to be an important prior form 
for the following reasons^ 


1) It is an improved representation of uncertainty over the 
2— point prior of the prior knowledge because it has 5 
adjustable parameters, namely the three possible values of 
the lot quality, and the prior probabilities associated with 
any two of them (the third is then known) - 

2) The analytical solution for the meeting point can be obtaine 
here as shown in the Section below. 
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3) The solution of the 3— point prior model enables us to draw 
some inferences about the behaviour of models based on 
k— point prior distributions with k 5: 4. 



lot quality (p) *■ 

FIGURE 3.2.2. The indifference lot quality, p^ . 

We define the indifference quality p^ (Figure 3.2.2) as the lot 
quality at which the two terminal decisions (accept or reject the 

lot) will be equally preferred (when the lot quality is p = p^, the 

risks in accepting or rejecting the lot are equal). If the three 
p^*s in the three point prior model are all on one side of the 
indifference lot quality p^, the decision to be made becomes 
trivial: one would then either always accept or always reject the 
lot regardless of its quality. Also, if no p^ exists such that 0 i 
p < 1, it only means that the losses W (p) and W .(p) have not been 

correctly quantified for a rational basis for accepting/rejecting 

the lot to exist. 

For the present purpose, we will impose the following restrictions 
on the selection of p^ , 


Pg and p^: 
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P-| < P 3 ^ Pg Pq Pe' (3.2.1) 

If we additionally locate at p^, that provides much computational 

convenience. In reality, we find that the restriction p_ = p can 

O o 

be relaxed only slightly because as we have found, there is a 
restriction on the magnitude of the difference between and p : if 
I p “ p I exceeds a certain value, the meeting point does not exist- 

w O 

For notational convenience, we also define the differential loss, 

W, (p), as U (p) -W (p). 
a a r 

3.2.1 THE MEETING POINT 

The Equations (2.2.7) and (2.2.8) reduce to the following form for 
the 3— point prior distribution 

- r n-r 

\ W,(p.) - c/p. p. a. (p.) (1 - p. ) "^ "^ = 0 (3. 2. 1.1) 

Z d ^1 ^1 1 *^1 '^1 

i=l'' 

r n — r 

S W,(p.) a. (p.) (1 - p. ) "* "* = 0 (3. 2. 1.2) 

Z d 1 

i=1 

A special case of this situation was solved by Wetherill C19I] who 
showed that the above equations assume a tractable form if one 
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assumes that the prior has a form such that p (1 - p ) ^ for k 

K 

=1^2 and 3- The equations then become solvable by the Horner^s 
method. However , we find this restriction to be unnecessary. One is 
able to obtain a closed form analytical solution for the meeting 
point as follows. 


We begin with the following notational 


q . ( 1 - p. ) 




kg 4— (a^/a^) 




-1 


-1 


“3 ‘V^i’ (“d'l’a’ "z ■ Pi 

kd <*3/a,) [Wd(P 3 > P3 - c][Md(p^) p, 

'1 - [ph«i / Pid2] [P3 / qj 


(Pa"! ^ Pi'^a] " [P 3 ' "'1] 


n 


m 


substitutions 




-1 


-1 


he above substitutions reduce equations (3.2. 1.1) and <3.2. 1.2) to 
wo linear equations in and Xg 

+ k^ X^ + kg Xg = 0 and (3.2. 1.3) 


^ k3 Xi . k4 Xg 


0 


(3.2. 1.4) 



III. 10 


resulting in the solution 


(k^ - 

•< 4 ) 

<k.k. - 
1 4 

'<e''3> ’’ 

(3.2. 1.5) 

"^3 - 

ki) 

<k.k- - 
1 4 

kgkg)--' 

(3.2. 1.6) 


Recall that -Ck^l are -functions o-f the loss -functions, the sampling 
cost and the prior distribution.) 


Back— substituting for and and replacing the above solutions 

(equations 3.2. 1.5/6) by k and k. respectively, 

A u 


[Pgq, / [P2 ^ qj ” = 

(03 ^ -<1) = 


(3.2. 1.7) 


(3.2. 1.8) 


Taking logarithms on both sides of equations (3.2. 1.7) and 
( 3 . 2 . 1 . 8 ), we obtain the following solution for the meeting point 

n^ = In k^ ~ ‘^b ^"[Pg^l'^Pl^ J]H (3. 2. 1.9) 

— -I 

"m “ “b ]H (3.2.1.10) 

where , 




In[p3qi/Piq3) ln[q2/qi] “ I"[p2q/Piq2] »"[P3^Pl] 


(3.2.1.11) 
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3.E.2 CONDITIONS FOR A MEETING POINT TO EXIST 

For a meeting point exist, equations (3.2. 1.9) and 

(3.2.1.10) must have a feasible solution. In addition, n and r 

m in 

must be greater than O. The necessary conditions for a meeting 
point to exist will now be discussed. 

Condition 1: A > 0 

It is easily seen from (3-2.1.11) that A is zero if and only if p^ = 
Pg “ P^- Thus a condition A = 0 in a real decision problem will 
never occur- 


Condition 2: k , k* > O 
a b 

This condition is required to be satisfied because both k and k. 

a b 

are arguments to logarithms in the equations for n^ and r^- From 


equations 

(3. 2. 1.5) 

and (3.2. 1.6) 

, it 

can be 

seen 

that 

(k^ 

k^ ) , 

(k3 - k^) 

and (k.k. 

1 4 

- k^k^) all 

must 

have 

the 

same 

sign 

The 

resulting 

conditions 

in terms of 

the 

losses , 

the 

cost 

of 

sampl ing 


and the prior distribution are easy to derive, and we will here give 
only an illustration of what these conditions imply. When p^ = 
the only constraint to be satisfied for Condition 2 to hold is that 




(3.2.2. 1 ) 
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The above inequality places a constraint on the "resolution" of the 

sampling plan: for a given differential loss function W^(p) and the 

cost of sampling c, the exact meeting point cannot be obtained if 

the difference between Pg and p^ is less than what is specified in 

the above equation. The larger the differential loss Wj(p) and the 

smaller the cost of sampling, the "finer" is the resolution of the 

sequential plan in terms of the possible values of p , p and p„ 

I ^ 3 • 


Condition 3: n , r >0 
m m 

We start with making a proposition [P about A which we will prove 
later in this section: 

Pr-opcnfittcfix IP : A is positive for p^ < p^ <Pg. 

Because we consider our problem only for the case where p^ < p^ < Pg 
(see Equation (3.2.1)), we are assured that the denominator in the 
equations for n and r is positive. 

Analytically finding the condition for the numerators to also be 
positive turns out to be a mathematically very complicated problem. 
We will discuss some of the numerical results obtained in this 
regard. Taking the cue from the results of Condition 2, for Pg = 
p^ , if we try to find a relationship between ( Pg — p^ ) and 
cri/W^(p^) + 1/W^(Pg)l, we find that 
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<Po - P^> > 1.5 C 


1/W^(p^) 


+ 1/W^(P2>! 


(3.S.2.2) 


guarantees that n_> O and r > 0, if both 

n> m 

greater than < Pg ” p^)/(1.5 c). For 

= <P 2 - P^)/(1.5 c), both 

equal to zero. 


W . ( p . ) and 
d ^1 

the special 

n and r 
m m 


Wd<P 2 ^ are 
case where 
are exactly 


Px-oof of Propositton tP: 

From equation (3.2.1.11), we -see that A = O at P-| ~ P 3 and Pg = Pg. 
It follows that A(p ) has at least one extremum for p^ < p„ < p_. On 

O I O 

partially differentiating A w.r.t p^, we also find that A has only 
one extremum for 0< P 2 < 1. Thus A has only on& extremum for p^ < p^ 

< Pg. Now, if and denote the slope of A at pg = p^ and Pg = p^ 

respectively, then 

= ^ln(1 - Pg) - ln(1 - P^i^jpg ~ P^ - iTi Pe}^^ " Pe^ (3. 2. 2. 3) 

^g = [ln(1 - Pg) - ln(1 - P^)]p^ - [in p^ - In Pgj ( 1 - p^) (3. 2. 2. 4) 

because there is only one extremum for A(pg) for p^ < Pg < Pg, A ^>0 
mplies Ag<0 and vice versa. From Equations (3. 2. 2. 3) and (3. 2. 2. 4), 
re get ^ O, as pg > p^. This implies that the slope of A at 

.g = p^ is positive and hence A itself is positive for P^ < P 3 p 2 - 



CHAPTER IV 


SENSITIVITY ANALYSIS 


4.0 INTRODUCTION 


Optimal (sequential) sampling plans (OSPs) obtained from the k— point 
prior distribution minimize the expected (Bayes) risk of sequential 
decisions in acceptance sampling as guided by the sequential 
sampling plans obtained using the theory given in Chapters II and 
III. In this chapter we will study the effect of changing the 
various cost, loss and prior distribution parameters on these 
sampling plans. In addition, we would also study the effect of 



knowledge of these parameters 

on 

the 

optimal i ty 

of 

the 

sampling 

plan, as shown in Figure 1.2. 

2 . 

Ne 

will 

be 

here 

be 

concerned 

only with the 3-poirU prior case. 

for 

whi ch 

the 

meet ing 

point can 

be obtained analytically. It 

is 

our belief 

that 

gene ral 


k— point prior distributions will give similar results. 

The following cases will be studied in this chapter: 

(1) The effect of the cost of inspection, on 

(a) The expected risk of using the OSP . 

(b) The ASN curve for the OSP. 

(c) The OC curve for the OSP. 
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t2) The effect of the choice of the prior distribution, on 

(a) The expected risk of using the OSP. 

(b) The DC curve of the OSP. 


(3) The effect of using a meeting point different from the 
exact meeting point, on 

(a) The decision boundaries of the OSP. 

<b) The expected risk of using the OSP. 


(4) The 
the 


effect of imperfect knowledge of the 
expected risk of using the resulting 


inspection cost, on 
sub-optimal plan. 


(5) The 
the 


effect of imperfect knowledge of the decision losses, on 
expected risk of using the resulting sub-optimal plan. 


(6) The 
the 


effect of imperfect inspection, on 
expected risk of using the resulting 


sub— optimal plan. 


4.1 THE MODEL 

While, in general, any arbitrary cost model can be used for 
generating OSPs for the k— point prior case, we will here assume a 
simple case where rejected items are recycled at a fixed unit cost. 
The monetary loss of accepting a lot is only due to the defective 
items that remain in it after inspection, the loss of accepting each 
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such (defective) item being the expected loss due to "external 
failure" (if this item is to be used in— house), assumed here to be a 
constant. The cost of inspecting one item is modeled as a constant. 


The decision losses and the inspection cost thus become: 


U (p) 
ac 

W . (p) 
rj ^ 


= the expected loss of accepting the lot 
= the expected loss of rejecting the lot 
= the cost of sampling 


L p C 


EF" 


•- Secyc' 


where 


^EF 

= the 

cost 

of 

external failure of one item. 

^RECYC 

= the 

cost 

of 

recycling one item. 

C 

s 

= the 

cost 

of 

inspecting one item, and. 

L 

= lot 

size . 




At the irvdi//<sr-&yxc<9' quality level, p^, the losses resulting from 
decisions "accept lot" and "reject lot" are equal and the decision 
maker is indifferent between choosing either of these two decisions. 
In the above model, the indifference quality is 


Po “ ^^RECYC ^ ^Ef] 


The "ideal" loss function is thus (from Equation 2.6.1): 
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4.3 THE IDEAL LOSS FUNCTION 


From Equation E.6.1, the ideal loss function for the above problem 
is s 


W .(p) 
rj ^ 



lot quality (fraction defective) p * 


FIGURE 4.1.1. The "ideal" loss function. 


An "ideal" loss function can only be realized when the lot quality 
can be known perfectly, without incurring any inspection cost. 


4.3 NUMERICAL RESULTS FOR AN OPTIMAL SSP . 
(Tabulated in TABLE A.1 in Appendix A.) 


In this section we present the numerical results obtained for an 
optimal SSP for the following example: 
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Losses and inspection cost: 

C^P = 100 

^RECYC 

C =1 

s 

thus, the indifference quality level = 10/100 = 0.1 


Prior 

Pi 

Ps 

P3 

^1 


, 1/3 
*2 

0.2 

P2 

FIGURE 4.3.1. The 3-point prior distribution. 


probabi 1 i ty 


1/3 


1/3 

a„ 


0.04 


0.1 


distribution (FIGURE 4.3.1): 

= 0.04 

= 0.2 

= 0.1 (fixed at p^, the indifference quality level) 

^2 ~ *3 ~ (i.e., p^, Pg and p^ are equally likely) 


Using the analytical results for the meeting point in Section 3.2.1, 
the meeting point of the acceptance and the rejection boundaries for 
this example is obtained at (202.590, 20.330). The effective 
sample size, N* (Section 2.3.1c), is subsequently obtained at a 
sample size of 176, where (176, 17) is an acceptance point and (176, 
18) a rejection point. The decision boundaries obtained using the 
algorithm BOUNDARIES are shown in Figure 4.3.2. 
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0 1 2 3 4 5 < 7 « 10 11 12 13 1*1 15 1 <- 17 

number of items sampled (n) 


FIGURE 4.3.2. 

sample size of 176, where (176, 17) is an acceptance point and (176, 
18) a rejection point. The decision boundaries obtained using the 
algorithm BOUNDARIES are shown in Figure 4.3.2. 

Once the above boundaries have been obtained, the calculation of the 
ASN and OC curves can proceed as stated in Section 2.5.2, using 
equations 2.5.2. 1 and 2. 5. 2. 2. The ASN and the OC curves are shown 
in Figures 4.3.3 and 4.3.4 respectively. 


The expected loss of using the OSP shown in Figure 4.3.2 can now be 
obtained using the results of Section 2.6. The expected loss 
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of using the OSP is compared with the "ideal" loss in Figure 4.3.5. 
The loss due to imperfect information, <5, which is also the loss of 
wrong decision, is shown in Figure 4.3.6. 

4.4 SENSITIVITY ANALYSIS 

C 1 ) Effect of the cost of inspection. 

(la) Expected Risk. (Figure 4.4.1.) 

The expected risk goes down as the cost of inspection 
decreases. Reduced cost of inspection permits the decision 
maker to obtain a larger sample size, and thus improve his 
knowledge of the lot quality. (TABLE A. 2 in Appendix A.) 

(lb) The ASN curve. (Figure 4.4.2.) 

The Average Sampling Number increases with decreasing cost of 
inspection because information about the lot quality can be 
obtained at a lesser cost. (TABLE A. 3 in Appendix A.) 

(lc) The OC curve. (Figure 4.4.3.) 

The DC curve shifts towards the ideal because more information 
about the lot quality is available at a lesser cost. 

(TABLE A. 3 in Appendix A.) 
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tflOO 

EFFECT OF THE PRfOR DISTRIBUTION QM . pa) AND INSFECTibN COST c ON RISK 
joo (p 3 IS ATTHE INDIFFERENCE OIMLITT LEVEL. 0 . 1 ) 
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FIGURE 4.4.1 
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FIGURE 4.4.E 
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LOT QUALITY (p) 


FIGURE 4.4.3. 


(2) Effect of the choice of the prior distribution. 


(2a) The expected risk. (Figure 4.4.1) 

At a lower value of the cost of inspection (c = 1)» the 
expected risk is more sensitive to the choice of the prior 
distribution than at higher values of the cost of inspection. 
(TABLE A. 2 in Appendix A.) 

(2b) The OC curve. (Figure 4.4.4) 

The OC curve becomes more discriminating as | p^ “ P.jl becomes 
smaller, as suggested by intuition. But at the same time, the 
average sampling number increases. (TABLE A. 3 in Appendix A.) 
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LOT QUALITY (p) 


FIGURE 4.4.4. 

3) The effect of using a different meeting point 

(Figures 4.4.5 and 4.4.6) j (TABLE A. 6(b) in Appendix A.) 

3a) The decision boundaries. (Figure 4.4.5.) 

If the exact meeting point cannot be calculated, optimal 
decision boundaries can still be obtained by assuming the 
meeting point to exist at a sufficiently large sample size. In 
fact, it is possible to obtain optimal decision boundaries if 
the meeting point is assumed to exist at a sample size that is 
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FORCeD MEETING POINT 



ON THE DECUaONBOUN?AFaE8 

-5 

0 5cy 400 » 20) 

NLf;©ER OF ITEMS 8A>^LED (n) 


FIGURE 4.4.5. 


greater than or equal to the effective sample size, N . If the 
meeting point is forced at (N^, R^), then should be chosen 
as a sufficiently large number, close to to facilitate 
numerical computations. The above result is a very important 
one in that it provides for an optimal solution for the OSP 
even if the meeting point cannot be calculated. 


If, on the other hand, the meeting point is forced at a sample 
size less than N*, the total oxpootod Bayes risk (averaged over 
the prior distribution) increases, even though the expected 
risk now is lesser in the vicinity of p^ (because of a smaller 
sample size). See Figure 4.4.6. (TABLE A. 6(a) in Appendix A.) 
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LOT QUALmr (py) 


FIGURE 4.4.6 

(4) The effect of imperfect knowledge of the inspection cost on the 
expected risk (Figure 4.4.7.): 

As B. result of imperfect informeition about the cost of 
inspection, the decision boundaries obtained will be different 
from the optimal decision boundaries. Again, the totcLt 

risk of the plan obtained using imperfect information will be 
higher, even though for a partial set of values of 
quality, the expected risk will turn out to be lesser. 

(TABLE A. 4 in Appendix A.) 
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FIGURE 4.4.11. 


(5) The effect of imperfect knowledge of the decision losses on the 
expected risk (Figures 4.4.8, 4.4.9 and 4.4.10.): 

(Comments same as in (4).) 

(TABLE A. 5 in Appendix A.) 


(6) 


The effect of imperfect inspection on the expected risk 
(Figures 4.4.11 and 4.4.12)- 


It is observed 
Type II errors, 
results in much 


from the comparison of the effects of Type I and 
that Type I error (rejection of a good item) 
larger risks than the Type II error. This is 



IV. 18 



FIGURE 4.4.12 

obvious because in a lot of very good quality, even a very /G\jt> 
incorrectly rejected items will lead to the rejection of the 
&rLtix'& lot. The obvious implication of this result is the need 
for accurate inspection equipment and well trained inspectors. 
In addition, an item should preferably be rejected only after a 
r e— ins pe ct ion , which will minimize the possibility of a Type I 
error. Also, in any appraisal of the "quality" of inspection, 
it should be well understood that the two types of error have 
vastly different monetary implications. 

(TABLE A. 7 in Appendix A.) 



CHAPTER V 


CONCLUSIONS AND DIRECTIONS FOR FUTURE WORK 


5.1 CONCLUSIONS 

A inocl 0 l W3S developed To obTdin opTintdl sequenTial seinvpling pl^ns 
(optimized on <e'xp<ec xish) based on (1) a generalized discrete 
prior distribution, called the k— point prior distribution, and (2) 
the concept of the "meeting point" of decision boundaries separating 
regions of different decision preferences. Various techniques for 
locating the meeting point so as to yield optimal plans were 
demonstrated: (1) numerical solution for k ^ 4, (2) analytical 

solution for k = 2 and k = 3, and (3) forcing a meeting point at a 
large enough sample size for all k ^ 2. The concept of the 

"expected risk of a sampling plan" was introduced and the effect of 
various input factors on the expected risk and the optimality of the 
plans were studied- 

We find that the strength of the k-point prior model lies in the 
it permits of modeling decision situations which can 
potentially have any arbitrary loss functions, using the most 
general (discrete) representation of prior information. Most 

Dx-ctc: t tl <2 clL problems can be solved satisfactorily using only the 
3-point prior distribution which was sol ved exotc t ly to obtain the 


meeting point. 
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It was also found 

(for the example 

considered) 

that 

the 

optimum 

(btained is a very 

“robust" optimum 

in the sense 

that 

the 

maximum 

idditional loss of 

under— (or over—) 

estimating 

the unit 

losses and 


lost-s by B,s ffiuch b.s 2!S/^ is only dboiit i% in th© wox^st ccl^& ■ On "th© 
3 ther hand, a s©©mingly “small" Type I inspection error of 5X might 
listurb the optimum by nearly 507, in the worst case. Inspection 
Brror of Type II offsets the effect due to the Type I error but only 
/ery marginally. Type II inspection errors are in comparison not 
very important as long as they are not too large. 

5.2 APPLICATION OF THE MODEL TO INSPECTION PERFORMANCE. 

Inspection error, notably of Type I (rejection of a non-defective 
item) can very adversely affect the optimality of the sampling plan, 
as is evident from sensitivity analysis in Chapter IV. We can apply 
the present (k—point prior) model to optimally track and improve 
inspection performance, targeting one type of inspection error at a 
time, as follows. 

Suppose we target the Type I error first. For this we set up an 
inspector audit scheme (i.e., we inspect the inspection process 
itself) in which the targeted inspector (or inspection station) is 
fed only good items, "n" representing the cumulative count these 
good items, without the inspector being aware of it. The number of 
such items called bad (or defective) by the inspector is the reject 
count "r" and this is plotted on the (n,r) sample space. 
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Assuming 

that the inspector is 

making Type I errors 

with a frequency 

of <p Tor 

every 

unit sar^pled, there are two different actions 

that 

could be 

taken : 

(1) correction 

of the error level 

from 4> — *■ 


where 

is the 

minimum error 

level possible; the 

additional 

cost 

involved 

here 

is the cost 

of retraining the 

inspector 

(or 


re— cal i brat i on of the inspection equipment), and, (2) maintenance of 
a status quo, i.e., no retraining is done and the error level 
remains at Suppose the loss of maintaining a Type I inspection 
error level of over a single day (or any other convenient time 
unit) of inspection is a function L(<^). Then, the loss caused by 
action (1) is 




(5.2.1 ) 


a constant, and the loss caused by action (2) is 


W2(<f!>) = L(<^) , 


(5.2.2) 


which will be usually be a monotoni cal ly increasing function of 4> . 


Additionally, if the cost of sequentially auditing one instance of 
inspection is known to be c, then the optimization problem becomes 
that of determining >j>h.sthe-r to retrain the inspector or not, 
minimizing the expected loss of the consequences of this decision. 
If the a priori knowledge of the unknown inspection error rate <p is 
expressed as a prior distribution, it is easily seen that the 
problem reduces to exactly the same form as described in Section 1.2 
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and can be solved by representing the prior distribution as a 
c-point prior distribution and using the methodology described in 
:hapter II. An optimal solution for sequentially tracking the Type 
LI error level can be developed on similar lines. 

>.3 DIRECTIONS FOR FUTURE WORK. 

’erhaps the most critical aspect of the construction of optimal 
sequential sampling plans is a reliable quantification of the 
lecision losses and the inspection cost. A method for 
kssessment of losses and costs accruing to the department or to the 
msiness unit (for which sampling inspection is carried out) should 
>e developed, perhaps as an accounting information system, which 
hould be sensitive to any changes within and outside the 

irganization which affects the losses/costs. Better methods for 
dnpointing the location of the (risk-) indifference quality level 
ire needed. 

inalytical solutions to the k— point prior model for k ^ 4 will be 

iseful in obtaining a greater insight into the problem, even though 
t is doubtful that these solutions will result in any real saving 
n terms of computational effort. This is because the method of 
orcing the meeting point at a large enough sample size will (we 
■resume) almost always be successful in obtaining optimal decision 
oundaries, though this has yet to be established. 
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Group sampling procedures have been completely overlooked in this 

I ^ ;t i <; clear that the optimal solution to these procedures 
work, ana it 

^ solving a very general scenario in acceptance 

will amount to =■ ►' 

, Ki,,tP that group sampling procedures will be better than 

sampling. Note 

.. , _.r>iina procedures only if the cost of sampling is a 

sequential sampH^y e f y 

» ♦ with respect to the number of items sampled in a 

concave function wi r h 

group (of 2: 1 items)* 
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APPENDIX A 


DATA TABLES FOR SENSITIVITY ANALYSIS 



TABLE A.1 ; RESULTS FOR THE EXAMPLE IN SECTION 4.3 TABLE A.2: EFFECT OF PRIOR DISTRIBUTION AND 

INSPECTION COST ON RISK 
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TABLE A-4; EFFECT OF IMPERFECT KWOWLEDGE 
OF THE INSPECTION COST, « 


1 


lot 

ADDITIONAL 

quality 

;; 

iP ) 

c = 1.0 

0.01 

' 0.00 

0.02 

, 0.00 

0.03 

0.00 

0.04 

0.00 

0.05 

0.00 

0.06 

0.00 

0.07 

0.00 

0.08 

0.00 

0.09 

0.00 

0.10 

0.00 

0.11 

0.00 

0.12 

0.00 

0.13 

0.00 

0.14 

0.00 

0.15 

0.00 

0.16 

0.00 

0.17 

0.00 

0.18 

0.00 

0.19 

0.00 

0.20 

0.00 

0.21 

0.00 

0.22 

0.00 

0.23 

0.00 

0.24 

0.00 

0.25 

0.00 


WHEN TRUE c - 1 


c = 0.5 

c « 2.0 

4.41 

- 4.44 

4.76 

- 4.40 

3.88 

- 1.03 

- 2.11 

11.05 

- 17.43 

33.30 

- 36.35 

53.83 

- 41.97 

54.20 

- 25.26 

30.44 

0.99 

- 0.05 

14.96 

- 14.93 

7.22 

- 5.77 

- 14.32 

19.22 

- 35.28 

45.26 

- 46.24 

61.63 

- 46.98 

66.34 

- 38.34 

59.15 

- 28.05 

47.81 

- 18.40 

36.39 

- 10.84 

24.35 

- 5.56 

16.65 

- 2.16 

9.33 

- 0.10 

4.99 

1.09 

213 

1.72 

C.31 

2.04 

- 0.82 



TAOi A.5: EFFECT OF IMPERFECT KNOWLEDGE OF LOSSES 
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TABLE A- 6 (a) 


t»ms 

(”) 


1 

2 

5 

4 

g 
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87 

88 
88 

70 

71 

72 

73 
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87 
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TABLE 4-6(b) 


lot 

quality 

(P) 


0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.10 

0.11 

0.12 

0.13 

0.14 

0.15 

0.16 

0.17 

0.18 

0.19 

0.20 

0.21 

0.22 

0.23 

0.24 

0.25 


EFFECT OF INEXACT LOCATION OF THE 
MEETING POINT 

ADDITIONAL RISK OF LOCATING THE 
MEETING POINT AT : 


(176.17) 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


( 148 , 14 ) 


0.00 

0.00 

0.03 

0.45 

2.24 

5.61 

8.30 

7.56 
3.66 
-0.48 
-2.40 
-1.91 
-0.35 

0.99 

1.56 
1.49 
1.11 
0.70 
0.39 
0.20 
0.09 
0.04 
0.02 
0.00 
0.00 


(119.11) 


-0.08 

- 0.20 

0.65 

7.41 

26.52 
51.87 
63.16 
48.08 

17.41 
-6.79 
-11.71 
-0.52 
15.43 
26.66 
29.90 
26.74 

20.52 
14.00 
8.67 
4.93 
2.58 
1.22 
0.49 
0.13 
-0.04 
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0.22 0.00 9.52 -10.62 -12.76 26.99 -15.24 -17.61 0.22 

0.23 0.00 6.49 -7.47 -9.32 18.10 -11.28 -13.52 0.23 

0.24 ‘ 0.00 4.63 -5.51 -7.15 12.56 -8.71 -10.84 0.24 

0.25 0.00 3.48 -4.24 -5.72 9.11 -6.97 -9.00 0.25 



PROGRAM LISTINGS 


Program BOUNDARIES 


: Program to find the (1) Effective Meeting Point, 
boundaries for the case of a 3 — point prior. 


and (E) The Decision 


‘ine MAXN 1000 

•ine MAXR 1000 

lefine maximum (computationally) possible size of the sample spaces a 
-ectangular integer grid <;(0,0), (0, MAXR), (MAXN,0), (MAXN, MAXR) J 


;lude <stdio.h> 

include standard I/O header file */ 


:lude <math.h> 

include math header file for calculating powO, abs ( ) */ 
it th1, thE, th3, a1, aE, a3, pol, poE, po3; 

:th1, thE, th3) are the 3 points on the prior. NOTE: thi > th3 > thE */ 
(a1, aS, a3) are the corresponding prior probabilities */ 

[pol, poS, po3) hold the posterior values for (a1, aE, a3) ■»/ 

at ral , rrl, raE, rrE, ra3, rr3, el, eE, e3; 

(ral, raE, ra3) are the losses (risks) of accepting lot at (thi, thE, 
)*/ 

(rrl, ..., ...) ... . . rejecting ... 

(el, eE, e3> defined in program #/ 
at ext, recyc, insp, ks; 

axt = cost of external failure, recyc = recycling cost, */ 

ks = inspection cost. ALL COSTS ARE per item. */ 


at RCMAXN3EMAXR3; 
risk matrix */ 


max_r , lotsize, 
max_r is the defect 
r = defect count, n 
array dib holds the 
dEb 


r, n, dIbCMAXNO, dEbCMAXNO; 
count at the meeting point, */ 

= number of items sampled */ 

defect count at all n for the REJECTION boundary.*/ 
ACCEPTANCE . */ 


E ♦out; 
ascii file 


contains meeting point, decesion boundaries */ 


ern double acc(); 
ern double rej ( ) ; 
ble ininrisk( ) ; 
ble posterior(); 
functions used */ 


i n ( > 


FILE ♦inp, ^extf; 
inp = f open( "inp", "r") ; 
out = f open ( "outf " , "w" ) ; 
extf = f open ( "inpg" , "w" ) ; 

fscanf(inp, '“/.d V.i 7.i /CfXnZf Xf \n'',&lotsize , &ext, Srecyc, &ks, Sa1 

2 ) ; 

insp *- lotsize; cost #= lotsize; ext *= lotsize; recyc *= lotsize; 
/■» costs pertaining to the entire lot */ 

th3 = ( recyc )/< ext ) j 

/ -» break-even or risk— indifference quality level */ 

pr intf ( "indif f erence quality is X— 5.3f \nEnter p2, p1", th3); 
scanf ( "■/.f/.f ■■,&th2, &th1); 

ra1 = acc(th1); ra2 = acc(thE); ra3 = acc(th3); 
rr1 = rej(th1); rr2 = rej(th2); rr3 = rej(th3); 

printf <" 7.-7. 3f 7-7. 3f 7.-7.3f \n7.-7.3f 7.-7. 3f 7.-7.3f \n7.-7.3f Xn" , 
ra1, raE, ra3, rr1, rr2, rr3, ks>; 
e1 = ra1 - rr1; 

eE = rrE — raE; 

e3 = rr3 - ra3; 

a3 = 1-a2-a1; 

printf ("el = 7.-5. 2f eE = 7.-5. EfXn", e1 , eE); 
meeting_point ( ) ; 

/* function finds meeting point */ 
boundaries ( ) ; 

/* function finds decision boundaries */ 

/* main */ 


ieting_point ( > 

* finds the meeting point for the 3— point prior */ 

float k1, kE, k3, k4, kden, ka, kb, aa, ab, ba, bb; 
float nr_den, real_n, real_r; 

k1 = - (a2/a1 )*(e2/e1 ) ; 
kE = (a3/a1 )*( e3/e1 ) ; 

l <;3 = - (a2/a1 )*(e2*th2 + ks)/(e1*th1- ks); 

k4 = (a3/a1 )*(e3*th3 - ks)/(e1*th1 - ks); 

kden = (k1*k4 - k2*k3); 
ka = (k3 - kD/kden; 
kb = (kE - k4)/kden; 


if 


((ka <= 0)!!(kb <= 0)) f ^ 

pr intf ( "Negative log— domain- ka-4 7.4f 

exitd) ; 


kb=7.-7.4f\n". 


/* if ka or kb <= 0, the logarithms cannot be evaluated */ 


ka, I 


aa = (th2/th1) * ( 1 -th 1 )/( 1 -thE) ; 
ah = (th3/th1) * (1-th1)/(1-th3) ; 


ba 

bb 


(1-th2)/(1-th1); 
(1-th3)/(1-th1) ; 


nr_den = ( log (ab )*log ( ba )-log (aa)^Hog ( bb) ) ; 

/* the denominator of the expression for real_n and real_r 

real_n = ( 1 og ( kb )-»log (ab)-log (aa)»log ( ka) ) /nr_den ; 
real_r = ( log ( ba >«log ( ka)-log ( bb)*log ( kb) ) /nr_den ; 
n = ceil(real_n) ; r = floor ( real_r ) ? 

/* the meeting point is converted into an integer point */ 

pr intf ( "n=X— 7 .3f r=X— 7.3f\n", real_n, real_r); 

/* meeting point */ 

/ 

nclude "risk.h" 


uble posterior (x,y) 

finds posterior at (n,r)=(x,y) given 3-point priors al, a2, a3 */ 
computationally the most expensive procedure */ 
t X, y ; 


double tot, r1, r2, r3; 

r1 = a1 * pow(th1 , (double)y) * pow( 1-th 1 , (double )( x-y) ) ; 

r2 = a2 * pow( th2, ( double )y) * powC 1-th2, (double) (x-y )) j 

r3 = a3 * pow ( th3, ( double )y ) * pow( 1—th3, (double) (x—y) ) ; 

tot = r1+r2+r3; 

po1 = rl/tot; 

po2 = r2/tot; 

po3 = 1— pol— po2; 

/* posterior */ 


uble minrisk(x, y, MIN, TAG) 

finds the decision with the minimum risk at point (x,y) and */ 
assigns the minimum risk value to MIN and the decision type to TA(5 


t 

t 

oat 


X, y ; 

*TAG; 

*MIN; 


tot, r1, r2, r3; 

R1, R2, R3j 

* pow( th 1 ,( double )y ) 
» pow( th2, (double )y) 

* pow ( th3 ,{ double )y ) 
r 1+r2+r3; 

r 1/tot ; 


double 
double 
r 1 = a1 

r2 = a2 

r3 = a3 

tot = 
pol = 

po2 = r2/tot; 
po3 = 1-po1-po2; 

/* calculate the values for the posterior 
R1 = po1*rr1 + po2*rr2 + po3«-rr3; 

R2 = po1*ra1 + po2'»ra2 + po3*ra3; 

R3 = ks + po1*(th1*RCx+1DCy+13+(1-th1)*RCx+13Cy3) 

+ po2-»(th2»RCx+13Cy+13+(1-th2>*RCx+1ICy3) 

+ po3*(th3*Ri:x+13Cy+13+(1-th3)*RCx+13Cy3); 

»MIN = (R1<R2)?( (R1<=R3)?(R1 ) = <R3> 1 < <R2<-R3)?(R2) ! (R3) ) 

*TAG = (R1<R2)?((R1<=R3)?(1)*(3)):((R2<=R3)?(2):(3>)j 

/* tag puts a "tag" on the decision with min risk */ 

/* roinrisk ♦/ 


pQw( 1-th1, (double) (x-y) ) j 
pow( 1-th2, (double) (x-y) ) j 
pow( 1-th3, (double ) (x-y) ) ; 


distribution at (x. 


V) 





indar i es ( ) 

char ch; 

int CNT, remain, cents, thous, i, j, k, max_n, max_y, d1, dS, t 

f 1 oat min , RSK ; 

d1 bCnD=r ; 
d2blln3=r+1 ; 

/* boundaries at the meeting point #/ 

for ( i=0; i<=r— 1 ;++i ) { 

posterior (n , i ) ; 

RCn3CiD = po1*ra1 + po2*ra2 + po3*ra3; 

> 

for (i = r; i<=n;++i)-C 

posterior (n , i ) ; 

RCnDCiD = po1*rr1 + po2»rr2 + po3*rr3; 

> 

/* the risks at meeting point— x (n,xx) are found #/ 


max_n = n ; 
max_r = r+1; 

/■X- meeting point is initially at n */ 


1 


/ * 


for (j=n-1jj>=0; — j)f 

d2 = -1; d1 = 0; 

for ( i=0; i<=max_r ;++i ) -C 

minrisk (j, i, &RSK, &tag ) ; 
RCjlCin = RSK; 

d2 = (tag == 2) ? (++d2) : (d2); 
d1 = (tag == 3) ? (++d1) : (d1); 

d2bCj3=d2; 
d1bCjD=d2+d1+1; 
max_r = d1bCj3; 

max n = (d1 == O) ? (j) s (max_n); 


printf ("\n") ; 

f printf (out, ‘•■/.-5d 7,-5d\n" ,max_n , d1bCmax_nd ) ; 
n = max_n; 

/* update meeting point •»/ 

pr intf ( '‘/C-3d 7 - 5 d\n",max_n ,d2bCmax_n3) ; 

-.- 5 ^ .-10.3^ 

i ,d2bEiD,d1bi;iD,d1bCi3-d2b[;i3-1, RCi3Cd2bCi3S , 


RCi3Cd1bi:i 


/# order: sample number; boundaries d2, d1; #of C points 

f close ( out ) ; 
boundaries */ 



2) Program RISK 




^ ASN calculation — output of mp used as input ^ 

^ OC curve calculation, as a spin off of ASN calculation ^ 

•5?^ RISK calculation, from the knowledge of OC curve and ASN ^ 

^ Data transfer between program mp and this one is through file **out" ^ 

— 


^define MAX 1000 

l^inc 1 ude <stdio«h> 


“ILE ^Qutr *»^inp, ^asnout; 

“out*® is the output of program mp, and is the input for this program 

"loat rrCMAXn, raCMAX3, ext, recyc, cost, insp; 
int d1bCMAX3, d2bCMAX3, n, r, lot; 

dlbCil is the rejection boundary at i */ 
dSbCil is the acceptance boundary at i */ 

(n,r) is the revised meeting point 


nain ( ) 

r 

int i , z j 

/* z is there just to absorb un-needed data in the input file "ou 
float p, incr, start, end; 

out=fopen( "outf ", "r" ) ; 
inp = f open( "inp", "r">; 
asnout = fopenC'af", “w"); 
fscanf (out,"7.d ■/.d\n",&n, 4r); 
for (i = 1 ; i<=n; ++!)•£ 

fscanf (out, "7.d 7.d 7.d Xd Xf Xf\n", &z, &d2bi:iD,&d1bCi D, 8,2 

SraCil, ArrCiD); 

dlbCiH += 1 ; d2bCi3 += 1 ; 

/* because the index of an array in C cannot be negative 

> 

fclose(out); . ♦ « \ 

fscanf (inp, "Xd Xf Xf Xf Xf, Slot, Sext, Srecyc, Scost, Sinsp) , 

054 ^ #= lot; recyc *= lot; cost *= lot; insp *= lot; 
printf( "start, end, incrementXn" ) ; 
scanf ("XfXfXf ", Sstart, Send, Sincr); 
for(p=start; p<=end; p+=incr ) ASN(p); 
f close(asnout) ; 

> 

ttinclude "risk.h" 


/* this function calculates the ASN at quality level p */ 
/« and as a by product, the OC curve too */ 
float p; 


t 


float total, cum_asn = 0 . 0 , 

float prCMAXlCMAXO, Pa=0.0, 

/* prCxOCyO = probability of 
sampling ♦/ 


probCMAXD, q, tmp; 
R, AC, RJ, MINRSK; 
reaching point (x,y) 


during sequent 


/» Pa is the probability of acceptance at quality level p */ 


int i T j ; 

q = '^-p; 

for (j=0; j<=10; ++ j ) prCODC j 3=0 . 0 ; 

pr CODE 1 3=1 . 0 ; ^ 

d1bC03=d1bC13; 

d2bC03=0 ; 

total = 0.0; 

Pa = 0.0; 

/* initialization over */ 

for <i = 1; i< =n ; ++ i ) { 
probEi3 = 0.0; 

/* initialize ■»/ 

prCi3C03 = 0; 

/# "boundary condition" «/ 

for (j=1; j<=d1bCi3; ++ j ) f 

prCi3Cj3 = q*prCi-13Cj3 + p^fprCi-13Cj-13 ; 

> 

/* iterative equation for REACHABLE boundary points •»/ 
j = d1b[Ii3; 

/* saves typing effort */ 

if ((j > i + 1) I I ( j > d1bCi-13)) prCi3Cj3 = 0.0; 

/* if rejection boundary unreachable */ 

for (j=0; j<=d2bCi3-1; ++ j ) prCi3Cj3 = 0.0; 

/* for unreachable points •»/ 

j = d2bCi3; 

if (j > 0)f 

if (d2bCi3 <= d2bCi-13) prCi3Cj3 = 0.0; ‘ 

/* if acceptance boundary unreachable */ 

/-» otherwise ... */ 

Pa += prCi3Cj3; 
probCi3 = prCi3Cj3; 
prCi3Cj3 = 0.0; 

/* j < O are unreachable points, stay clear from them! */ 

j = d1bCi3; 
probCi3 += prCi3U3; 
prCi3Cj3 = 0.0; 

/* probC3 accumulates probabilities at ALL termination 
points for the current lot quality p ♦/ 

total += probCi3; 
cum_asn += i* pr'obCi3; 

> /* repeat for all sampling points */ 



RJ = r ej ( p) ; 

MINRSK = (AC <= RJ) ? (AC) : (RJ); 

R = Pa * AC + ( 1 - Pa) * RJ ; 

printfC'p = 7.-7. 3f total = 7.-7. 4f ASN = 7.-7. Of Pa = 7.-7. 5f 
R = 7— 10.3f\n‘', p, total, cum_asn/total , Pa/total, R) ; 
f printf (asnout, "7.-7. 3f 7.-7. Of 7.-7. 5f 7.-10. 3f 7.-10. 3f 7.-10. 3f\n", 
cunri_asn. Pa/total, MINRSK, R-MINRSK, R+cuiri_asn*iii5p/lot-MINRSK ) 



DISCUSSION 


The following points were raised during the thesis examinations 

(1) At the theoretical meeting point, one is indifferent between 
either accepting or rejecting the lot, or taking an additional 
observation and subsequently making a terminal acce pt /re j ect 
decision. In real life, should one r-arLdorril’^ select a terminal 
decision or should one take one more sample? 

Commj&rt t s : 

In reality, the decision boundaries will not meet because the 
meeting point (n , r ) will in general be real (See Figure 4.4.5.). 
So the situation of randomly selecting a terminal decision will not 
arise. If such a situation arises, one would probably take one more 
sample and end up conclusively in one of the two terminal decision 
regions . 


(2) Can a generalization be made that inspection errors will otlwctys 
contribute heavily to sub~opt imali ty of the sampling plan? 


CarrmuBTits: 


No, such a generalization cannot be made. But for a few 
representative problems we considered, inspection errors did 
contribute heavily to sub-optimality. Whether this happens or not 
depends on the loss/cost structure of the problem. 



curve more 


(3) Will a decreasing j - p^| always make the OC 
discr i mi Dating? 

Camnu&rtt B : 

A decrease in j p^ — p^j may not always discriminate between p^ and 
pg better, but will always discriminate between "good" and "bad" 
lots better. 

(4) What are the possible values of the posterior distribution at 
the meeting point, when the starting prior for the 3-point case is 

1/3 1/3 1/3 

0.04 6.1 0.2 p > 

Comrru&rxt £ : 

If the risk-indifference lot quality is at 0.1, the posterior 
probability of p = 0.1 will be high as compared to other values of p 
at the meeting point: 


0.04 0.1 0.2 

(5) How does the "additional information" curve for lot quality p 
vary with the number of items sampled? 

CorwfnjBTLts: 

As sampling proceeds, the information about p decreases. 

Thus the above curve will be a monotoni cal ly decreasing curve. 




